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Relatively uniform light curves and spectral evolution of Type la supernovae (SNe la) have led 
to the use of SNe la as a "standard candle" to determine cosmological parameters, such as the 
Hubble constant, the density parameter, and the cosmological constant. Whether a statistically 
^ ■ significant value of the cosmological constant can be obtained depends on whether the peak 

\^ ' luminosities of SNe la are sufficiently free from the effects of cosmic and galactic evolutions. 

OO . Here we first review the single degenerate scenario for the Chandrasekhar mass white dwarf 

CO ' (WD) models of SNe la. We identify the progenitor's evolution and population with two chan- 

C ; nels: (1) the WD-hRG (red-giant) and (2) the WD-hMS (near main-sequence He-rich star) 

■ channels. In these channels, the strong wind from accreting white dwarfs plays a key role, which 
yields important age and metallicity effects on the evolution. 

We then address the questions whether the nature of SNe la depends systematically on envi- 

■ ronmental properties such as metallicity and age of the progenitor system and whether significant 
Oh' evolutionary effects exist. We suggest that the variation of the carbon mass fraction X{C) in 

the C-l-0 WD (or the variation of the initial WD mass) causes the diversity of the brightness 
5^ ' of SNe la. This model can explain the observed dependence of SNe la brighness on the galaxy 

^ ! types. 

. Finally, applying the metallicity effect on the evolution of SN la progenitors, we make a 

prediction of the cosmic supernova rate history as a composite of the supernova rates in different 
types of galaxies. 

1. Introduction 

Type la supernovae (SNe la) are good distance indicators, and provide a promising tool 
for determining cosmological parameters (e.g., Branch 1998). SNe la have been discov- 
ered up to z ~ 1.32 (Gilliland et al. 1999). Both the Supernova C osmology Project (Perl- 



mutter et al. 1997, 1999) and the High-z Supernova Search Team ( Garnavich et al. 1998 ; 
Riess et al. 1998) have suggested a statistically significant value for the cosmological 
constant. 

However, SNe la are not perfect standard candles, but show some intrinsic variations 
in brightness. When determining the absolute peak luminosity of high-redshift SNe la, 
therefore, these analyses have taken advantage of the empirical relation existing between 
the peak brightness and the light curve shape (LCS) . Since this relation has been obtained 
from nearby SNe la only (Phillips 1993; Hamuy et al. 1995; Riess et al. 1995), it is 
important to examine whether it depends systematically on environmental properties 
such as metallicity and age of the progenitor system. 

High-redshift supernovae present us very useful information, not only to determine 
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cosmological parameters but also to put constraints on the star formation history in the 
universe. They have given the SN la rate at z ~ 0.5 (Pain 1999) but wiU provide the SN 
la rate history over < z < 1. With the Next Generation Space Telescope, both SNe la 
and SNe II will be observed through z ~ 4. It is useful to provide a prediction of cosmic 
supernova rates to constrain the age and metallicity effects of the SN la progenitors. 

SNe la have been widely believed to be a thermonuclear explosion of a mass-accreting 
white dwarf (WD) (e.g., Nomoto et al. 1997a for a review). However, the immedi- 
ate progenitor binary systems have not been clearly identified yet (Branch et al. 1995). 
In order to address the above questions regarding the nature of high-redshift SNe la, 
we need to identify the progenitors systems and examine the "evolutionary" effects (or 
environmental effects) on those systems. 

In §2, we summarize the progenitors' evolution where the strong wind from accreting 
WDs plays a key role (Hachisu, Kato, & Nomoto 1996, 1999, hereafter HKN96, HKN99; 
Hachisu, Kato, Nomoto & Umeda 1999, HKNU99). In §3, we addresses the issue of 
whether a difference in the environmental properties is at the basis of the observed range 
of peak brightness (Umeda et al. 1999b). In §4, we make a prediction of the cosmic 
supernova rate history as a composite of the different types of galaxies (Kobayashi et al. 
1999). 



2. Progenitor's evolution and the white dwarf wind 

There exist two models proposed as progenitors of SNe la: 1) the Chandrasekhar mass 
model, in which a mass-accreting carbon-oxygen (C-t-0) WD grows in mass up to the 
critical mass Mja ~ 1.37 — 1.38M0 near the Chandrasekhar mass and explodes as an SN 
la (e.g., Nomoto et al. 1984, 1994), and 2) the sub-Chandrasekhar mass model, in which 
an accreted layer of helium atop a C-l-0 WD ignites off-center for a WD mass well below 
the Chandrasekhar mass (e.g., Arnett 1996). The early time spectra of the majority 
of SNe la are in excellent agreement with the synthetic spectra of the Chandrasekhar 
mass models, while the spectra of the sub-Chandrasekhar mass models are too blue to 
be comparable with the observations (Hoflich & Khokhlov 1996; Nugent et al. 1997). 

For the evolution of accreting WDs toward the Chandrasekhar mass, two scenarios have 
been proposed: 1) a double degenerate (DD) scenario, i.e., merging of double C+0 WDs 



with a combined mass surpassing the Chandrasekhar mass limit (Iben & Tutukov 1984 



Webbink 1984), and 2) a single degenerate (SD) scenario, i.e., accretion of hydrogen-rich 
matter via mass transfer from a binary companion (e.g., Nomoto 1982; Nomoto et al. 
1994). The issue of DD vs. SD is stiU debated (e.g., [Branch et al. 1995), although the- 



oretical modeling has indicated that the merging of WDs leads to the accretion-induced 
collapse rather than SN la explosion (Saio & Nomoto 1985, 1998; segretain et al. 1997). 

In the SD Chandrasekhar mass model for SNe la, a WD explodes as a SN la only 
when its rate of the mass accretion (M) is in a certain narrow range (e.g., Nomoto 1982; 
Nomoto & Kondo 1991). In particular, if M exceeds the critical rate Mb in Eq.(2.1) 
below, the accreted matter extends to form a common envelope (Nomoto et al. 1979). 
This difficulty has been overcome by the WD wind model (see below). For the actual 
binary systems which grow the WD mass (Mwd) to Mja, the following two systems are 
appropriate. One is a system consisting of a mass-accreting WD and a lobe-filling, more 
massive, slightly evolved main-sequence (MS) or sub-giant star (hereafter "WD-I-MS 
system"). The other system consists of a WD and a lobe-filling, less massive, red-giant 
(hereafter "WD-I-RG system"). 
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Figure 1. Ratio of the photospheric velocity to the escape velocity there Uph/uesc is plotted 
against the decreasing rate of the envelope mass for WDs with masses of O.GMq, Q.7Mq, O.SMq, 
0.9Mq, I.OMq, 1.2Mq, 1.3Mq, and 1.377Mq. We regard the wind as "strong" when the 
photospheric velocity exceeds the escape velocity there, i.e., Vph > Vesc- If not, it is regarded as 
"weak." 



2.1. White dwarf winds 

Optically thick WD winds are driven when the accretion rate M exceeds the critical rate 
Mb. Here Afb is the rate at which steady burning can process the accreted hydrogen into 
helium as 

Mb « 0.75 X 10-^ - 0.40^ Mq yr-\ (2.1) 

With such a rapid accretion, the WD envelope expands to i?ph ~ O.IRq and the 
photospheric temperature decreases below logTph 5.5. Around this temperature, the 
shoulder of the strong peak of OPAL Fe opacity ( Iglesias fc Rogers 1993[ ) drives the 
radiation-driven wind (HKN96; HKN99). We plot the ratio I'ph/i'csc between the pho- 
tospheric velocity and the escape velocity at the photosphere in Figure |l| against the 
mass transfer rate. We call the wind strong when Vph > Vcsc- When the wind is strong, 
Vpii ~ 1000 km s~^ being much faster than the orbital velocity. 

If the wind is sufficiently strong, the WD can avoid the formation of a common envelope 
and steady hydrogen burning increases its mass continuously at a rate Mb by blowing 
the extra mass away in a wind. When the mass transfer rate decreases below this critical 
value, optically thick winds stop. If the mass transfer rate further decreases below ~ 0.5 
Mb, hydrogen shell burning becomes unstable to trigger very weak shell flashes but still 
burns a large fraction of accreted hydrogen. 

The steady hydrogen shell burning converts hydrogen into helium atop the C-l-0 core 
and increases the mass of the helium layer gradually. When its mass reaches a certain 
value, weak helium shell flashes occur. Then a part of the envelope mass is blown off but 
a large fraction of He can be burned to C-l-0 (Kato & Hachisu 1999) to increase the WD 
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Figure 2. An illustration of the WD-j-RG (symbiotic) channel to Type la supernovae. 



mass. In this way, Thus strong winds from the accreting WD play a key role to increase 
the WD mass to Mia. 

2.2. WD+RG system 

This is a symbiotic binary system consisting of a white dwarf (WD) and a low mass 
red-giant (RG). A full evolutionary path of the WD+RG system from the zero age main- 
sequence stage {stage A) to the SN la explosion {stage F) is as follows. 
(A) Zero age main-sequence. 
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(B) The primary has evolved first to become an asymptotic giant branch (AGB) star 
and blows a slow wind (or a super wind) at the end of the AGB evolution. 

(C) If the superwind from the AGB star is as fast as or slower than the orbital ve- 
locity, the wind outflowing from the system takes away the orbital angular momentum 
effectively. As a result the wide binary shrinks greatly (by about a factor of ten or more) 
to become a close binary. This is a similar process to the common envelope evolution. 

(D) Then the AGB star undergoes a common envelope evolution. The AGB star forms 
a C+O WD, and the initial secondary remains a main-sequence star (MS). 

(E) The initial secondary evolves to a red-giant (RG) forming a helium core and fills 
up its inner critical Roche lobe. Mass transfer begins. The WD component blows a 
strong wind and the winds can stabilize the mass transfer even if the RG component has 
a deep convective envelope. 

(F) The WD component has grown in mass to Mia and explodes as a Type la super- 
nova. 

The occurrence frequency of SNe la through this channel is much larger than the earlier 
scenario, because of the following two evolutionary processes, which have not considered 
before. 

(1) Because of the AGB wind at stage C, the WD -I- RG close binary can form from 
a wide binary even with such a large initial separation as < 40, OOOi?©. Our earlier 
estimate (HKN96) is constrained by < l,5OOi?0. 

(2) When the RG fills its inner critical Roche lobe, the WD undergoes rapid mass 
accretion and blows a strong optically thick wind. Our earlier analysis has shown that 
the mass transfer is stabilized by this wind only when the mass ratio of RG/WD is smaller 
than 1.15. Our new finding is that the WD wind can strip mass from the RG envelope, 
which could be efficient enough to stabilize the mass transfer even if the RG/WD mass 
ratio exceeds 1.15. If this mass-stripping effect is strong enough, though its efficiency rycff 
is subject to uncertainties, the symbiotic channel can produce SNe la for a much (ten 
times or more) wider range of the binary parameters than our earlier estimation. 

With the above two new effects (1) and (2), the WD-fRG (symbiotic) channel can ac- 
count for the inferred rate of SNe la in our Galaxy. The immediate progenitor binaries in 
this symbiotic channel to SNe la may be observed as symbiotic stars, luminous supersoft 
X-ray sources, or recurrent novae like T CrB or RS Oph, depending on the wind status. 

2.3. WD+MS system 

In this scenario, a C-l-0 WD is originated, not from an AGB star with a C-l-0 core, 
but from a red-giant star with a helium core of ^ 0.8 — 2.0Mq. The helium star, which 
is formed after the first common envelope evolution, evolves to form a C-l-0 WD of 
~ 0.8 — I.IMq with transferring a part of the helium envelope onto the secondary main- 
sequence star. 

As an example for this system, let us consider a pair of 7A/q -t- 2.5Mq with the initial 
separation of ~ 50 — 6OOi?0. The binary evolves to SN la through the following stages 
(stages A-F in Figure | and G in HKNU99): 

(1) stage A-C: When the mass of the helium core grows to I.OMq < Mi^ylc < 1.4Mq, 
the primary fills its Roche lobe and the binary undergoes a common envelope evolution. 

(2) stage C-D: After the common envelope evolution, the system consists of a helium 
star and a main-sequence star with a relatively compact separation of a/ ^3 — 40 Rq 
and Forb 0.4 - 20 d. 

(3) stage D: The helium star contracts and ignites central helium burning to become a 
helium main-sequence star. The primary stays at the helium main-sequence for ^ 1 x 10^ 

yr- 
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Figure 3. An illustration of the WD+MS channel to Type la supernovae through the 
common envelope evolution to the helium matter transfer. 



(4) stage E: After helium exhaustion, a carbon-oxygen core develops. When the core 
mass reaches 0.9 — I.OMq, the helium star evolves to a red-giant and fills again its inner 
critical Roche lobe. Almost pure helium is transferred to the secondary because the mass 
transfer is stable for the mass ratio q = M1/M2 < 0.79. The mass transfer rate is as 
large as ~ 1 x 10~^Mq yr~^ but the mass-receiving main-sequence star (~ 2 — 3Mq) 
does not expand for such a low rate. 

(5) stage F: The secondary has received 0.1 — O.4M0 (almost) pure helium and, as 
a result, it becomes a helium-rich star as observed in the recurrent nova U Sco. The 
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Figure 4. The region to produce SNe la in the log Po — Md,o plane for five initial white 
dwarf masses of O.75M0, 0.8Mq, 0.9Mq, I.OMq (heavy solid line), and I.IMq. The region of 
-^WD,o = 0.7Mq almost vanishes for both the WD+MS and WD+RG systems, and the region 
of MwD,o = O.75M0 vanishes for the WD+RG system. Here, we assume the stripping efficiency 
of rjcfi ~ 1. For comparison, we show only the region of A/wd.o ~ 1-OMq for a much lower 
efficiency of 77cff ~ 0.3 by a dash-dotted line. 

primary becomes a C+0 WD. The separation and thus the orbital period gradually 
increases during the mass transfer phase. The final orbital period becomes Porb ~ 0.5 — 40 
d. 

(6) stage G: The white dwarf accretes hydrogen-rich, helium-enhanced matter from a 
lobe-filling, slightly evolved companion at a critical rate and blows excess matter in the 
wind. The white dwarf grows in mass to Afja and explodes as an SN la. 

This evolutionary path provides a much wider channel to SNe la than previous sce- 
narios. A part of the progenitor systems are identified as the luminous supersoft X-ray 
sources (van den Heuvel et al. 1992) during steady H-burning (but without wind to avoid 
extinction), or the recurrent novae like U Sco if H-burning is weakly unstable. Actually 
these objects are characterized by the accretion of helium-rich matter. 

2.4. Realization Frequency 

For an immediate progenitor system WD+RG of SNe la, we consider a close binary 
initially consisting of a C+0 WD with Mwd.o = 0.6 — 1.2Mq and a low-mass red-giant 
star with A/rg,o = 0.7 — S.OAf© having a helium core of A^hcO = 0.2 — 0.46AfQ (stage 
E). The initial state of these immediate progenitors is specified by three parameters, i.e., 
AfwD.Oi -^^RG.o = -^d,Oi a-iid the initial orbital period Pq (Affjo.o is determined if Pq is 
given) . 

We follow binary evolutions of these systems and obtain the parameter range(s) which 
can produce an SN la. In Figure the region enclosed by the thin solid line produces SNe 
la for several cases of the initial WD mass, AfwD,o = 0.75 - 1.1 Mq. For smaller AfwD,0) 
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Figure 5. The regions of SNe la is plotted in the initial orbital period vs. the initial 
companion mass diagram for the initial WD mass of Mwd.o = I.OMq. The dashed and solid 
lines represent the cases of solar abundance {Z — 0.02) and much lower metallicity of Z = 0.004, 
respectively. The left and the right regions correspond to the WD+MS and the WD+RG 
systems, respectively. 



the wind is weaker, so that the SN la region is smaller. The regions of AfwD.o = O.GMq 
and O.TMq vanish for both the WD+MS and WD+RG systems. 

In the outside of this region, the outcome of the evolution at the end of the calculations 
is not an SN la but one of the followings: 

(i) Formation of a common envelope for too large Md or Pq ~ day, where the mass 
transfer is unstable at the beginning of mass transfer. 

(ii) Novae or strong hydrogen shell flash for too small Md.o, where the mass transfer 
rate becomes below 10~^ Mq yr~^. 

(iii) Helium core flash of the red giant component for too long Pq, where a central 
helium core flash ignites, i.e., the helium core mass of the red-giant reaches 0.46Mq. 

(iv) Accretion-induced collapse for MyvD.o > l.SM©, where the central density of the 
WD reaches ~ 10^" g cm~^ before heating wave from the hydrogen burning layer reaches 
the center. As a result, the WD undergoes collapse due to electron capture without 
exploding as an SN la (Nomoto & Kondo 1991). 

It is clear that the new region of the WD+RG system is not limited by the condition 
of q < 1.15, thus being ten times or more wider than the region of HKN96's model 
(depending on the the stripping efficiency of r]cs). 

The WD+MS progenitor system can also be specified by three initial parameters: 
the initial C+0 WD mass Mwd.Oj the mass donor's initial mass M(j,o, and the orbital 
period Pq- For Mwd,o = 1-OMq, the region producing an SN la is bounded by Md,o = 
1.8 — 3.2M0 and Pq = 0.5 — 5 d as shown by the solid line in Figure |^. The upper 
and lower bounds are respectively determined by the common envelope formation (i) 
and nova-like explosions (ii) as above. The left and right bounds are determined by the 
minimum and maximum radii during the main sequence of the donor star (IIKNU99). 
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Figure 6. WD mass vs. metallicity diagram showing the metalhcity dependence of optically 
thick winds We regard the wind as "strong" if the wind velocity at the photosphere exceeds 
the escape velocity but "weak" if the wind velocity is lower than the escape velocity. The 
term of "He" or "Fe" wind denotes that the wind is accelerated by the peak of iron lines near 
log T(K) ~ 5.2 or of helium lines near log T(K) ~ 4.6. The dashed line indicates the demarcation 
between the "strong" wind and the "weak" wind. 



We estimate the rate of SNe la originating from these channels in our Galaxy by using 
equation (1) of Iben & Tutukov (1984). The realization frequencies of SNe la through 
the WD+RG and WD+MS channels are estimated as - 0.0017 yr^^ (WD+RG) and 
- 0.001 yr-i (WD+MS), respectively. The total SN la rate of the WD+MS/WD+RG 
systems becomes ~ 0.003 yr~^, which is close enough to the inferred rate of our Galaxy. 

2.5. Low metallicity inhibition of type la supernovae 

In the above SN la progenitor model, the accreting WD blows a strong wind to reach 
the Chandrasekhar mass limit. If the iron abundance of the progenitors is as low as 
[Fe/H]< —1, then the wind is too weak for SNe la to occur. 

In this model, an interesting metallicity effect has been found (Kobayashi et al. 1998; 
Hachisu & Kato 1999). The wind velocity is higher for larger Mwd and larger Fe/H 
because of higher luminosity and larger opacity, respectively. In order to blow sufficiently 
strong wind {v^j > fesc), AfwD should exceed a certain mass (Fig. |l|). As seen from 
the dashed line in Figure ^, is larger for lower metalhcity; e.g., = 0.65, 0.85, 
and 0.95 Mq for Z = 0.02, 0.01, and 0.004, respectively. In order for a WD to grow its 
mass at M > Mb, its initial mass Mwd,o should exceed Mw. In other words, Mw is the 
metallicity-dependent minimum AfwD,o required for a WD to become an SN la. 

2.6. Possible detection of hydrogen 

In our scenario, the WD winds form a circumstellar envelope around the binary systems 
prior to the explosion, which may emit X-rays, radio, and Ha lines by shock heating 
when the ejecta collide with the circumstellar envelope. The mass accretion rate in the 
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present models is still as high as 1 x 10^® Mq yr~^ for some of the white dwarfs near 
the Chandrasekhar limit, so that such a white dwarf explode in the strong wind phase. 

Our strong wind model of case PI predicts the presence of circumstellar matter around 
the exploding white dwarf. Whether such a circumstellar matter is observable depends 
on its density. The wind mass loss rate from the white dwarf near the Chandrasekhar 
limit is as high as M '--^ 1 x 10^'* —1 x 10^^ A/q yr^^ and the wind velocity is w = 
1000 km s~^. Despite the relatively high mass loss rate, the circumstellar density is not 
so high because of the high wind velocity. For steady wind, the density is expressed by 
AI/v {— Anr'^p). Normalized by the typical red-giant wind velocity of 10 km s^^, the 
density measure of our white dwarf wind is given as AI/viq ^ 1 x 10^^" — Ix 10~^ Mq 
yr~^, where vio = v/10 km s~^. 

Behind the red-giant, matter stripped from the red-giant component forms a much 
dense circumstellar tail. Its rate is as large as ~ 1 x lO^^M© yr"^ with the velocity of 

100 km s~^. The density measure of the dense red-giant wind thus formed is M /vio 
1 X 10-^ Mq yr~i. 

Further out, circumstellar matter is produced from the wind from the red-giant com- 
panion, which is too far away to cause significant circumstellar interaction. 

For some cases, winds from the WD have stopped before the explosion. Therefore, 
circumstellar matter is dominated by the wind from the red-giant companion whose 
velocity is as low as ^ 10 km s^^. 

At SN la explosion, ejecta would collide with the circumstellar matter, which pro- 
duces shock waves propagating both outward and inward. At the shock front, particle 
accelerations take place to cause radio emissions. Hot plasmas in the shocked materials 
emit thermal X-rays. The circumstellar matter ahead of the shock is ionized by X-rays 
and produce recombination Ha emissions ( Gumming et al. 199(]| ). Such an interaction 
has been observed in Type lb, Ic, and II supcrnovae, most typically in SN 1993J (e.g., 
Suzuki & Nomoto 1995 and references therein). 

For SNe la, several attempts have been made to detect the above signature of cir- 
cumstellar matter. There has been no radio and X-ray detections so far. The up- 
per li mit set by X-ray obser vations of SN 1992A is M/vw = (2-3) x 10"^ Mq 
yr~^ ( Schlegel fc Petre 1993 ). Radio observations of SN 1986G have provided the most 
strict upper limit to the circumstellar density as M/vio = 1 x 10"'' Mq yr~^ (Eck et 
al. 1995). This is still 10 — 100 times higher than the density predicted for the white 
dwarf wind. If the WD wind has ceased and the wind mass loss rate from the red-giant 
is significantly higher than 1 x 10~^ Mq yi:~^, radio detection could be possible for 
very nearby SNe la as close as SN 1986G. (Note also that SN 1986G is not a typical SN 
la but a subluminous SN la.) 

For Ha emissions, Branch et al. (1983) noted a small, narrow emission feature at the 
rest wavelength of Ha, which is blueshifted by 1800 km s"^ from the local interstellar 
Ca II absorption. Though this feature was not observed 5 days later, such high velocity 
hydrogen is expected from the white dwarf wind model. For SN 1994D, Gumming et al. 
(1996) obtained the upper limit of M/vio = 6 x 10^^ Mq yr^^. Further attempts to 
detect Ha emissions are highly encouraged. 



3. The origin of diversity of SNe la and environmental effects 

There are some observational indications that SNe la are affected by their environ- 
ment. The most luminous SNe la seem to occur only in spiral galaxies, while both spiral 
and elliptical galaxies are hosts for dimmer SNe la. Thus the mean peak brightness is 
dimmer in ellipticals than in spiral galaxies (Hamuy et al. 1996). The SNe la rate per 
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Figure 7. Abundances in mass fraction in the inner core of the 6 Mq star for Y = 0.2775 
and Z = 0.02 at the end of the second dredge-up. 

unit luminosity at the present epoch is almost twice as high in spirals as in ellipticals 

(CappcUaro ct al. 1997). Moreover, Wang ct al. (1997) and Ricss ct al. (1999) found 
that the variation of the peak brightness for SNe located in the outer regions in galaxies 
is smaller. 

Hoflich et al. (1998, 1999) examined how the initial composition of the WD (mctallicity 
and the C/0 ratio) affects the observed properties of SNe la. Umeda et al. (1999a) 
obtained the C/0 ratio as a function of the main-sequence mass and metallicity of the 
WD progenitors. Umeda ct al. (1999b) suggested that the variation of the C/0 ratio 
is the main cause of the variation of SNe la brightness, with larger C/0 ratio yielding 
brighter SNe la. We will show that the C/0 ratio depends indeed on environmental 
properties, such as the metallicity and age of the companion of the WD, and that our 
model can explain most of the observational trends discussed above. We then make some 
predictions about the brightness of SN la at higher redshift. 

3.1. C'/O ratio in WD progenitors 

In this section we discuss how the C/0 ratio in the WD depends on the metallicity 
and age of the binary system. The C/0 ratio in C+0 WDs depends primarily on the 
main-sequence mass of the WD progenitor and on metallicity. 

We calculated the evolution of intermediate- mass (3 — 9Mq) stars for metallicity 
Z^O.OOl - 0.03. In the ranges of stellar masses and Z considered in this paper, the 
most important metallicity effect is that the radiative opacity is smaller for lower Z. 
Therefore, a star with lower Z is brighter, thus having a shorter lifetime than a star with 
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Figure 8. The total ^^C mass fraction included in the convective core of mass, M = 1.14Mq, 
just before the SN la explosion as a function of the C+O core mass before the onset of mass 
accretion, Mco ■ The lower bounds of Mco obtained from the age effects and the conditions for 
strong wind to blow are also shown by arrows. 

the same mass but higher Z. In this sense, the effect of reducing metallicity for these 
stars is almost equivalent to increasing a stellar mass. 

For stars with larger masses and/or smaller Z, the luminosity is higher at the same 
evolutionary phase. With a higher nuclear energy generation rate, these stars have larger 
convective cores during H and He burning, thus forming larger He and C-0 cores. 

As seen in Figure 0, the central part of these stars is oxygen-rich. The C/0 ratio is 
nearly constant in the innermost region, which was a convective core during He burning. 
Outside this homogeneous region, where the C-0 layer grows due to He shell burning, 
the C/0 ratio increases up to C/0 > 1; thus the oxygen-rich core is surrounded by a 
shell with C/0 > 1. In fact this is a generic feature in all models we calculated. The 
C/0 ratio in the shell is C/0 ~ 1 for the star as massive as ~ 7Mq, and C/0 > 1 for 
less massive stars. 

When a progenitor reaches the critical mass for the SNe la explosion, the central core 
is convective up to around 1.1 Mq. Hence the relevant C/0 ratio is between the central 
value before convective mixing and the total C/0 of the whole WD. Using the results 
from the C6 model (Nomoto et al. 1984), we assume that the convective region is 1.14 
Mq and for simphcity, C/0 = 1 outside the C-0 core at the end of second dredge- up. 
Then we obtain the C/0 ratio of the inner part of the SNe la progenitors (Fig. ^). 

From this figure we find three interesting trends. First, while the central C/0 is a 
complicated function of stellar mass (Umeda et al. 1999a), as shown here the C/0 ratio 
in the core before SNe la explosion is a decreasing monotonic function of mass. The 
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central C/0 ratio at the end of second dredge- up decreases with mass for Mms > 5Mq, 
while the ratio increases with mass for Mms ^4:Mq; however, the convective core mass 
during He burning is smaller for a less massive star, and the C/0 ratio during shell He 
burning is larger for smaller C+O core. Hence, when the C/0 ratio is averaged over 
1.1 Mq the C/0 ratio decreases with mass. Second, as shown in Umeda et al. 1999, 
although the C/0 ratio is a complicated function of metallicity and mass, the metallicity 
dependence is remarkably converged when the ratio is seen as a function of the C+O 
core mass (Mqq) instead of the initial main sequence mass. 

According to the evolutionary calculations for 3—9 Mq stars by Umeda et al. (1999a), 
the C/0 ratio and its distribution are determined in the following evolutionary stages of 
the close binary. 

(1) At the end of central He burning in the 3—9 Mq primary star, C/0< 1 in the 
convective core. The mass of the core is larger for more massive stars. 

(2) After central He exhaustion, the outer C+O layer grows via He shell burning, where 
C/0> 1 (Umeda et al. 1999a). 

(3a) If the primary star becomes a red giant (case C evolution; e.g. van den Heuvel 
1994), it then undergoes the second dredge-up, forming a thin He layer, and enters the 
AGB phase. The C+O core mass, Mqq, at this phase is larger for more massive stars. 
For a larger Mqo the total carbon mass fraction is smaller. 

(3b) When it enters the AGB phase, the star greatly expands and is assumed here to 
undergo Roche lobe overflow (or a super-wind phase) and to form a C+O WD. Thus the 
initial mass of the WD, Mwd,o, in the close binary at the beginning of mass accretion is 
approximately equal to Mco- 

(4a) If the primary star becomes a He star (case BB evolution) , the second dredge-up 
in (3a) corresponds to the expansion of the He envelope. 

(4b) The ensuing Roche lobe overflow again leads to a white dwarf of mass Mwd.o = 
Mco- 

(5) After the onset of mass accretion, the WD mass grows through steady H burning 
and weak He shell flashes, as described in the WD wind model. The composition of the 
growing C+O layer is assumed to be C/0=1. 

(6) The WD grows in mass and ignites carbon when its mass reaches Mia = 1.367M0, 
as in the model C6 of Nomoto et al. (1984). Because of strong electron-degeneracy, 
carbon burning is unstable and grows into a deflagration for a central temperature of 
8 X 10® K and a central density of 1.47 x 10^ g cm~'^. At this stage, the convective 
core extends to Mr = 1.14Mq and the material is mixed almost uniformly, as in the C6 
model. 

In Figure ||, we show the carbon mass fraction X{G) in the convective core of this 
pre-explosive WD, as a function of metallicity (Z) and initial mass of the WD before the 
onset of mass accretion, Mqo- Figure || reveals that: 1) X{C) is smaller for larger Mco — 
MwDfi. 2) The dependence of ^(C) on metallicity is small when plotted against Mco: 
even though the relation between Mco and the initial stellar mass depends sensitively 
on Z (Umeda et al. 1999a). 

3.2. Brightness of SNe la and the C/0 ratio 

In the Chandrasekhar mass models for SNe la, the brightness of SNe la is determined 
mainly by the mass of ^^Ni synthesized (MNise)- Observational data suggest that Mniss 
for most SNe la lies in the range MNise ~ 0.4 - 0.8Mq (e.g. Mazzali et al. 1998). This 
range of AfNise can result from differences in the C/0 ratio in the progenitor WD as 
follows. 

In the deflagration model, a larger C/0 ratio leads to the production of more nuclear 
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Figure 9. SN la frequency for a galaxy of mass 2 x 10 M© as a function of Mco for Z=0.02. 
For the WD+RG system, constraints from the companion's age are shown by the arrows. SNe 
la from the WD+MS system occur in spirals but not in ellipticals because of the age effect. 
Mco and MNise is assumed to be related as shown here. 



energy and buoyancy force, thus leading to a faster propagation. The faster propagation 
of the convective deflagration wave results in a larger MNise- For example, a variation of 
the propagation speed by 15% in the W6 - W8 models results in MNise values ranging 
between 0.5 and O.TMq (Nomoto et al. 1984), which could explain the observations. 

In the delayed detonation model, AfNise is predominantly determined by the deflagration- 
to-detonation-transition (DDT) density pddt, at which the initially subsonic deflagra- 
tion turns into a supersonic detonation (Khokhlov 1991). As discussed in Umeda et al. 
(1999b), pddt could be very sensitive to ^(C), and a larger ^(C) is likely to result in 
a larger pddt and MNise- 

Here we postulate that Mnisg and consequently brightness of a SN la increase as the 
progenitors' C/0 ratio increases (and thus Mwd,o decreases). As illustrated in Figure ^, 
the range of Mnisg ^ 0.5 — CSMq is the result of an X{C) range 0.35 — 0.5, which is the 
range of ^(C) values of our progenitor models. The A'(C) - Mnisg - Mwd,o relation we 
adopt is still only a working hypothesis, which needs to be proved from studies of the 
turbulent flame during explosion (e.g., Niemeyer & Hillebrandt 1995). 



3.3. Metallicity and age effects 
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3.3.1. Metallicity effects on the minimum MwD,o 

As mentioned in §2.5, is the metallicity-dependent minimum Mwd,o for a WD to 
become an SN la {strong wind condition in Fig. The upper bound Mwd,o — I.OTMq 
is imposed by the condition that carbon should not ignite and is almost independent of 
metallicity. As shown in Figure ^ the range of Mqq ~ Mwd.o can be converted into a 
range of ^(C). From this we find the following metallicity dependence for ^(C): 

(1) The upper bound of X(C), which is determined by the lower limit on Mqo imposed 
by the metalhcity-dependent conditions for a strong wind, e.g., X(C) < 0.51, 0.46 and 
0.41, for Z=0.02, 0.01, and 0.004, respectively. 

(2) On the other hand, the lower bound, ^(C) ~ 0.35 — 0.33, does not depend much 
on Z, since it is imposed by the maximum Mco- 

(3) Assuming the relation between MNise and X{C) given in Figure ||, our model 
predicts the absence of brighter SNe la in lower metallicity environment. 

3.3.2. Age effects on the minimum Mwd.o 

In our model, the age of the progenitor system also constrains the range of ^(C) in 
SNe la. In the SD scenario, the lifetime of the binary system is essentially the main- 
sequence lifetime of the companion star, which depends on its initial mass M2- HKNU99 
and HKN99 have obtained a constraint on M2 by calculating the evolution of accreting 
WDs for a set of initial masses of the WD (AfwD,o — Mqo) and of the companion (M2), 
and the initial binary period {Pq). In order for the WD mass to reach Mia, the donor 
star should transfer enough material at the appropriate accretion rates. The donors 
of successful cases are divided into two categories: one is composed of slightly evolved 
main-sequence stars with M2 ^ 1.7 — 3.6Af0 (for Z=0.02), and the other of red-giant 
stars with M2 ~ 0.8 - 3.1Mq (for Z=0.02) (Fig. |). 

If the progenitor system is older than 2 Gyr, it should be a system with a donor star 
of M2 < 1.7Mq in the red-giant branch. Systems with A/2 > 1.7Mq become SNe la in a 
time shorter than 2 Gyr. Likewise, for a given age of the progenitor system, M2 must be 
smaller than a limiting mass. This constraint on Af2 can be translated into the presence 
of a minimum Mqo for a given age, as follows: For a smaller Af2, i.e. for the older 
system, the total mass which can be transferred from the donor to the WD is smaller. 
In order for Af-yvD to reach Afja, therefore, the initial mass of the WD, AfvvD,o — AfcO) 
should be larger. This implies that the older system should have larger minimum Afco 
as indicated in Figure |[ Using the X(C)-Afco and Alj^i^g-X^C) relations (Fig. ||), we 
conclude that WDs in older progenitor systems have a smaller X{C), and thus produce 
dimmer SNe la. 

3.4. Comparison with observations 

The first observational indication which can be compared with our model is the possible 
dependence of the SN brightness on the morphology of the host galaxies. Hamuy et al. 
(1996) found that the most luminous SNe la occur in spiral galaxies, while both spiral 
and elliptical galaxies are hosts to dimmer SNe la. Hence, the mean peak brightness is 
lower in elliptical than in spiral galaxies. 

In our model, this property is simply understood as the effect of the different age of 
the companion. In spiral galaxies, star formation occurs continuously up to the present 
time. Hence, both WD-I-MS and WD-I-RG systems can produce SNe la. In elliptical 
galaxies, on the other hand, star formation has long ended, typically more than 10 Gyr 
ago. Hence, WD-I-MS systems can no longer produce SNe la. In Figure ^, we show 
the frequency of the expected SN I for a galaxy of mass 2 x lO-'^^Af© for WD-I-MS and 
WD-I-RG systems separately as a function of Afco- Here we use the results of HKN99 



16 Nomoto, Umeda, Hachisu, Kato, Kobayashi, Tsujimoto: SN la 

and HKNU99, and the Mco -X{C) and Mnisg - ^(C) relations given in Figure |. Since 
a WD with smaUer Mco is assumed to produce a brighter SN la (larger MNise), our 
model predicts that dimmer SNe la occur both in spirals and in ellipticals, while brighter 
ones occur only in spirals. The mean brightness is smaller for ellipticals and the total 
SN la rate per unit luminosity is larger in spirals than in ellipticals. These properties 
are consistent with observations. 

The second observational suggestion is the radial distribution of SNe la in galaxies. 
Wang et al. (1997) and Riess et al. (1998) found that the variation of the peak 
brightness for SNe la located in the outer regions in galaxies is smaller. This behavior 
can be understood as the effect of metallicity. As shown in Figure || even when the 
progenitor age is the same, the minimum Mco is larger for a smaller metallicity because 
of the metallicity dependence of the WD winds. Therefore, our model predicts that the 
maximum brightness of SNe la decreases as metallicity decreases. Since the outer regions 
of galaxies are thought to have lower metallicities than the inner regions (Zaritsky et al. 
1994; Kobayashi & Arimoto 1999), our model is consistent with observations. Wang et 
al. (1997) also claimed that SNe la may be deficient in the bulges of spiral galaxies. This 
can be explained by the age effect, because the bulge consists of old population stars. 

3.5. Diversity of high redshift supernovae 

We have suggested that ^(C) is the quantity very likely to cause the diversity in Mwise 
and thus in the brightness of SNe la. We have then shown that our model predicts that 
the brightness of SNe la depends on the environment, in a way which is qualitatively con- 
sistent with the observations. Further studies of the propagation of the turbulent flame 
and the DDT are necessary in order to actually prove that X{C) is the key parameter. 

Our model predicts that when the progenitors belong to an old population, or to a 
low metal environment, the number of very bright SNe la is small, so that the variation 
in brightness is also smaller. In spiral galaxies, the metallicity is significantly smaller at 
redshifts z > 1, and thus both the mean brightness of SNe la and its range tend to be 
smaller. At z > 2 SNe la would not occur in spirals at all because the metallicity is too 
low. In elliptical galaxies, on the other hand, the metallicity at redshifts z ~ 1 — 3 is not 
very different from the present value. However, the age of the galaxies at 2; ~ 1 is only 
about 5 Gyr, so that the mean brightness of SNe la and its range tend to be larger at 
z > 1 than in the present ellipticals because of the age effect. 

We note that the variation of ^(C) is larger in metal-rich nearby spirals than in high 
redshift galaxies. Therefore, if X{C) is the main parameter responsible for the diversity 
of SNe la, and if the light curve shape (LCS) method is confirmed by the nearby SNe 
la data, the LCS method can also be used to determine the absolute magnitude of high 
redshift SNe la. 

3.6. Possible evolutionary effects 

In the above subsections, we consider the metallicity effects only on the C/0 ratio; this 
is just to shift the main-sequence mass - Mwd,o relation, thus resulting in no important 
evolutionary effect. However, some other metallicity effects could give rise to evolution 
of SNe la between high and low redshifts (i.e., between low and high metallicities). 

Here we point out just one possible metallicity effect on the carbon ignition density in 
the accreting WD. The ignition density is determined by the competition between the 
compressional heating due to accretion and the neutrino cooling. The neutrino emission 
is enhanced by the local Urea shell process of, e.g., ^^Ne-^^F pair (Paczyhski 1973). (Note 
that this is different from the coiivective Urea neutrino process). For higher metallicity, 
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dotted line is for the DD scenario, and the other lines are for our SD scenario with (solid 
line) and wit hout ("dashed line) the in etallicitv effect on SNe la. Observational data so urces: 



Nissen et al. (1994) ; open squares, Cratton (1991 



Hlled circles, Rdvard sson et al. (1993) ; open triang les. Barbuy & Erdclyi-Mcndcs (1989) 



stars, 



the abundance of ^^Ne is larger so that the cooling is larger. This could delay the carbon 
ignition until a higher central density is reached (Nomoto et al. 1997d). 

Since the WD with a higher central density has a larger binding energy, the kinetic 
energy of SNe la tends to be smaller if the same amount of ^^Ni is produced. This might 
cause a systematically slower light curve evolution at higher metallicity environment. 
The carbon ignition process including these metallicity effects as well as the convective 
Urea neutrino process need to be studied (see also Iwamoto et al. 1999 for nucleosynthesis 
constraints on the ignition density). 



4. The chemical evolution in the solar neighborhood 

The role of SNe II and SNe la in the chemical evolution of galaxies can be seen in 
the [O/Fe]- [Fe/H] relation (Fig. |l^: Metal-poor stars with [Fe/H] < -1 have [O/Fe] 
~ 0.45 on the average, while disk stars with [Fe/H] > — 1 show a decrease in [O/Fe] with 
increasing metallicity. To explain such an evolutionary change in [O/Fe] against [Fe/H], 
we use the chemical evolution model that allows the infall of material from outside the 
disk region. The infall rate, the SFR, and the initial mass function (IMF) are given by 
Kobayashi et al. (1998), and the nucleosynthesis yields of SNe la and 11 are taken from 
Nomoto et al. (1997bc) and Tsujimoto et al. (1995). 

The metallicity effects on SNe la on the chemical evolution is examined. For the 
DD scenario, the distribution function of the lifetime of SNe la by Tutukov & Yungelson 
(1994) is adopted, majority of which is ~ 0.1 — 0.3 Gyr. Figure |l^ shows the evolutionary 
change in [O/Fe] for three SN la models. The dotted line is for the DD scenario. The other 
lines are for our SD scenario with (solid line) and without (dashed line) the metallicity 
effect on SNe la. 
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(i) In the DD scenario the hfethiie of the majority of SNe la is shorter than 0.3 Gyr. 
Then the decrease in [0/Fe] starts at [Fe/H] ^ ~2, which is too early compared with the 
observed decrease in [0/Fe] starting at [Fe/H] ^ — 1. 

(ii) For the SD scenario with no metaUicity effect, the companion star with M ~ 
2.QMq evolves off the main-sequence to give rise to SNe la at the age of ~ 0.6 Gyr. The 
resultant decrease in [0/Fe] starts too early to be compatible with the observations. 

(iii) For the metaUicity dependent SD scenario, SNe la occur at [Fe/H] > — 1, which 
naturally reproduce the observed break in [0/Fe] at [Fc/H] ^ —1. Also the low- 
metallicity inhibition of SN la provides a new interpr etation of the SN Il-lik e abundance 
patterns of the Galactic halo and the DLA systems ( [Kobayashi et al. 1998 ). 



5. Cosmic supernova rates 

Attempts have been made to predict the cosmic supernova rates as a function of redshift 



by using the observed cosmic star formation rate (SFR) (Ruiz-Lapuente fc Canal 1998 



3adat et al. 1998; Yimgelson fc Livio 1998). The observed cosmic SFR shows a peak 



at z ~ 1.4 and a sharp decrease to the present (Madau et al. 1996). However, UV 
luminosities which is converted to the SFRs may be affected by the dust extinction 



(Pettini et al. 1998). Recent updates of the cosmic SFR suggest that a peak lies around 
z - 3. 

Kobayashi et al. (1998) predicts that the cosmic SN la rate drops at z ~ 1 — 2, due to 
the metallicity-dependence of the SN la rate. Their finding that the occurrence of SNe la 
depends on the metaUicity of the progenitor systems implies that the SN la rate strongly 
depends on the history of the star formation and metal-enrichment. The universe is 
composed of different morphological types of galaxies and therefore the cosmic SFR is a 
sum of the SFRs for different types of galaxies. As each morphological type has a unique 
star formation history, we should decompose the cosmic SFR into the SFR belonging to 
each type of galaxy and calculate the SN la rate for each type of galaxy. 

Here we first construct the detailed evolution models for different type of galaxies 
which are compatible with the stringent observational constraints, and apply them to 
reproduce the cosmic SFR for two different environments, e.g., the cluster and the field. 
Secondly with the self-consistent galaxy models, we calculate the SN rate history for each 
type of galaxy and predict the cosmic supernova rates as a function of redshift. We adopt 
Hq = 50 km s~^ Mpc~^, fio = 0.2, Ao = and the galactic age of 15 Gyr for a standard 
model, which corresponds to the redshift at the formation epoch of galaxies Zf ~ 5. 

5.1. Supernova rates and galaxy types 

We assume that elliptical galaxies are formed by a single star burst and stop the star 
formation at t ^ 1 Gyr due to the supernova-driven galactic wind (e.g., Kodama & 
Arimoto 1997), while spiral galaxies are formed by a relatively continuous star formation. 
These models can well reproduce the present gas fractions, colors, supernova rates, and 
the color evolution in cluster ellipticals (see Figures |ll|and^. 

Present supernova rates observed in the various type of galaxies (Cappellaro et al. 
1997) put the constraints on the SN la progenitor models. Using the galaxy model 
shown in Figures |l| and m, we show that our SN la model well reproduces the present 
supernova rates in both spirals and ellipticals (Kobayashi et al. 1999). 

5.1.1. Spiral galaxies 

The observed SN II rate TZu in late-type spirals is about twice the rate in early-type 
spirals. On the other hand, the observed SN la rate Kia, in both types of spirals are 
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Figure 11. The star formation rates (top panel) gas fractions (middle panel) and the B — V col- 
ors (bottom panel) in four types of spirals : SOa-Sa (solid line), Sab-Sb (long-dashes line), Sbc-Sc 
(short-dashed hne), and Scd-Sd (dotted hne). The present B-V colors are taken from Roberts & 
Haynes (1994), and the present gas fractions are normal ized by the present blue luminosities an d 
are derived from HI fractions (Roberts fc Haynes 1994|) and H2/HI ratios (Casoli et al. 1998). 



nearly the same. Therefore the present TZii^/TZu ratio in early-type spirals is about twice 
that in late-type spirals. Such a difference in the relative frequency is a result of the 
difference in the SFR, because the dependences of TZu and T^ia on the SFR are different 
due to the different lifetimes of supernova progenitors. Therefore the observed 7?.ia/7^ii 
ratio gives a constraint on the SN la progenitor model. 

In the DD scenario, the lifetime of majority of SNe la is ^ 0.1 — 0.3 Gyr, so that the 
evolution of TZia is similar to that of TZu. Therefore TZia,/TZu is insensitive to the SFR. 
This results in the small differences in T^ia/T^ii among the various type of spirals, which 
is not consistent with observations. 

In our SN la model, if the iron abundance of progenitors is [Fe/II]> —1, the occurrence 
of SNe la is determined from the lifetime of the companions. If SNe la would occur only in 
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Figure 12. The passive color evolution of cluster ellipticals up to z < 
data are taken from Stanford et al. (1998). 



1. The observational 



the MS+WD systems with relatively short lifetimes, TZia,/TZu would have been insensitive 
to the SFR. On the contrary, if SNe la would occur only in the RG+WD systems with 
long lifetimes, the present difference in 7?,ia/7^ii between early and late type spirals would 
have been too large, reflecting the large difference in SFR at an early epoch. Owing to 
the presence of these two types of the progenitor systems in our SN la progenitor model, 
the observed difference in TZia/TZu can be reproduced. 

5.1.2. Elliptical galaxies 

We assume that a bulk of stars in cluster ellipticals are formed at z > 3, and have ages 
> 10 Gyr. In the DD scenario, the SN la lifetimes are mostly so short as tia < 0.3 Gyr 
that too few SNe la occur at the present epoch. Our SN la model includes the WD+RG 
systems with tia. > 10 Gyr, thus well reproducing the present SNe la rate in ellipticals. 

From z ~ 0.2 to z ~ 0, the SN la rate gradually decreases to the present by ~ 20%. 
This is due to the smallest companions with 0.9Mq to produce SNe la, which give the 
longest SN la hfetime of - 11 Gyr for Z = 0.002 and - 18 Gyr for Z = 0.02. As the star 
formation in ellipticals has been stopped more than 10 Gyr before, majority of SN la 
progenitors with Z = 0.002 have already explode and only metal-rich SNe la occurs at z ~ 
0. The decrease in the SN la rate from z 0.2 to z also depends on the star formation 
history in ellipticals and the galactic age. If ellipticals have undergone the relatively 
continuous star formation, as suggested by the hierarchical clustering simulations, the 
SN la rate might keep constant to the present. 

5.2. Cosmic supernova rates 

We calculate the cosmic supernova rate by summing up the rates of spirals and ellipticals 
with the ratio of the relative mass contribution among the types. The relative mass 
contribution is given by the observed relative luminosity proportion for ellipticals, SOa- 
Sa, Sab-Sb, Sbc-Sc, and Scd-Sd (Pence 1976) and the calculated mass to light ratio in 
B-band. We adopt the the initial comoving density of gas Ogoo — 3.5 x lO^^h"^. 
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Figure 13. The upper panel shows the cosmic SFR (sohd hne) as a composite of spi 
rals (l ong-dashed hnel an d ellipticals (sh ort-dashed line). The dots ar e the observational 



data ( |Galle go ct al. 1995|, open tria ngle; [Lilly et al. 1996| open circles; [Madau et al. 



open squares: 



1996 



Jpnnolly et al. 1997, filled circles: [Tres"se fc Maddox 199^ . open pentagon 



Treyer et al. 1995, filled square; CTazebrook et al. 19U^ , filled pentagon; Hughes et al. 1998 
star; Pettini et al. 1998, filled triangle). I'he lower panel shows the cosmic supernova rate (solid 
line) as a composite of spirals (long-dashed line) and ellipticals (short-dashed line). The obser- 
vational data are from Pain et al. (1996) and Pain (1999). 



5.2.1. In clusters 

First, we make a prediction of the cosmic supernova rates in the cluster galaxies using 
the galaxy models which are in good agreements with the observational constraints (Figs. 
|ll|and 12). The upper panel of Figure shows the cosmic SFR (solid line) as a composite 
of spirals (long-dashed line) and ellipticals (short-dashed line). 
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In our galaxy models, ellipticals undergo a star burst at z > 3 and the duration of the 
star formation is ^ 1 Gyr, while spirals undergo relatively continuous star formation. 
Thus, only the SFR in spirals is responsible for the cosmic SFR at z < 2. For reference, 
the observed cosmic SFR in field, so-called Madau's plot, is also plotted. Compared with 
them, the predicted cosmic SFR has a little shallower slope from the present to the peak 



at 2 1.4 (see also Totani et al. 1997). The high SFR in ellipticals appears at z > 3. 
Such ellipticals may be hidden by the dust extinction ( Pcttini et al. 199^ ) or may have 



formed at z 5 (Totani et al. 1997) 



The chemical enrichment in ellipticals has taken place with so much shorter timescale 
that the iron abundance reaches [Fc/H] ^ — 1 at z > 4 and the metallicity effect on SNe 
la does not appear. In spirals, the iron abundance reaches [Fe/H] ~ — 1 at z ~ 2.3. 

The lower panel of Figure |l^ shows the cosmic supernova rates (solid line) as a compos- 
ite of spirals (long-dashed line) and ellipticals (short-dashed line). The upper and lower 
three lines show the SN II and la rates, respectively. Our SN la model can successfully 



reproduce the observed rate at z ~ 0.5 ( Pain et al. 199€ , filled circle; Pain 1999 , filled 
square) . 

The SN la rate in spirals drops at z 2 because of the low-metallicity inhibition of 
SNe la. In ellipticals, the chemical enrichment takes place so early that the metallicity 
is large enough to produce SNe la at z > 3. Thus, the SN la rate depends almost only 
on the lifetime. A burst of SNe la occurs after ~ 0.5 Gyr from the beginning of the star 
formation which corresponds to the shortest lifetime of the WD-I-MS binaries; this forms 
a peak of the SN la rate at z 2.5. The second peak of the SN la rate appears at z ~ 1.5 
(~ 2 Gyr) due to the beginning of the explosions of SNe la from WD-I-RG binaries. If 
SNe la at z > 2 are observed with their host galaxies using the Next Generation Space 
Telescope, we can precisely test the metallicity effect by finding the drop of SN la rate 
in spirals. 

The cosmic SN la rate using the observed cosmic SFR drops at z ~ 1.6 (z 1.2 in 



Kobayashi et al. 1998). This is because the chemical enrichment of the whole universe 



is much slower than in spirals and ellipticals. 
5.2.2. Infields 

The observed spectra of ellipticals in the Hubble Deep Field suggest that the formation 
of ellipticals is protracted in fields, that is, the formation epochs of ellipticals span in 



the wide range of redshifts of 1 < z < 4 ( Pranceschini et al. 1998| ). We also predict 



the cosmic supernova rates for this case, assuming the distribution function of formation 
epochs to meet Franceschini's et al. (1998) results. 

The upper panel of Figure |lj shows the cosmic SFR (solid line) as a composite of 
spirals (long-dashed line) and field ellipticals (short-dashed line). 

The SFR in spirals is the same as in the upper panel of Figure ^ but the star formation 
in ellipticals continues to the present on the average. The peak of the star formation is 



at z ~ 3, which is consistent with the recent sub-mm data (Hughes et al. 1998). The 
synthesized cosmic SFR can successfully reproduce the observed one, except for the recent 
Ha data ( plazebrook et al. 1999 ). 



The lower panel of Figure |14| shows the cosmic supernova rates (solid line) as a com- 
posite of spirals (long-dashed line) and ellipticals (short-dashed line). The upper and 
lower three lines show the SN II and la rates, respectively. 

The SN la rate in spirals drops at z ~ 2 as in the lower panel of Figure |l3|. Contrary 
to Figure |l^, the SN la rate in field ellipticals gradually decreases from z ~ 2 to z 3, 
which is not due to the metallicity effect but due to the lifetime effect. Although the 
chemical enrichment timescale in each elliptical is as fast as in the lower panel of Figure 
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Figure 14. The same as Figure hSl but for the formation epochs of ellipticals spanning over 
1 < z < 4, which might correspond to field ellipticals. 



|l3|, the star formation averaged over whole ellipticals takes place more gradually, and the 
number of stars becomes smaller toward higher redshifts. 

In this case, SNe II may be observed in low-redshift eUipticals. The different SN II and 
la rates between cluster and field ellipticals reflects the difference in the star formation 
histories in different environments. 



5.2.3. Summary 

We make a prediction of the cosmic supernova rate history as a composite of the differ- 
ent types of galaxies. We adopt the SN la progenitor scenario including the metallicity 
effect proposed by Kobayashi et al. (1998), which can successfully reproduce the chem- 
ical evolution of the solar neighborhood. To calculate the comic SFR, we construct the 
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galaxy evolution models for spirals and ellipticals to meet the observational constraints 
such as the present gas fractions and colors for spirals, and the mean stellar metallicity 
and the color evolution from the present to z ~ 1 for ellipticals. 

Owing to the two types of the progenitor system (MS+WD and RG+WD), i.e., shorter 
(0.6 — 1 Gyr) and longer lifetimes (2 — 15 Gyr) of SN la progenitors, we can explain the 
moderate contrast of the relative ratio of SN la to SN II rate TZis,/TZu between the early 
and late types of spirals. Owing to over 10 Gyr lifetime of the RG+WD systems, SNe 
la can be seen even at the present in ellipticals where the star formation has already 
stopped over 10 Gyr before. 

Then we construct the cosmic SFR as the composite of the SFR for different types of 
galaxies, and predict the cosmic supernova rates: 

(a) In the cluster environment, the synthesized cosmic SFR has an excess at z > 3 
corresponding to the SFR in ellipticals and a shallower slope from the present to the peak 
at 2; ~ 1.4, compared with Madau's plot. The predicted cosmic supernova rate suggests 
that SNe la can be observed even at high redshifts because the chemical enrichment 
takes place so early that the metallicity is large enough to produce SNe la at 2 > 3 in 
cluster ellipticals. In spirals the SN la rate drops at z ~ 2 because of the low-metallicity 
inhibition of SNe la. 

(6) In the field environment, ellipticals are assumed to form at such a wide range of 
redshifts as 1 < z < 4. The synthesized cosmic SFR has a broad peak around 2: ~ 3, 
which well reproduces the observed one. The SN la rate is expected to be significantly 
low at 2 > 2 because the SN la rate drops at z ~ 2 in spirals and gradually decreases 
from 2: ~ 2 in ellipticals. 

This work has been supported in part by the grant-in- Aid for COE Scientific Research 
(07CE2002) of the Ministry of Education, Science, Culture and Sports in Japan. 
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